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Abstract 

The aim of the study was to compare the epidemic spread on static 
and dynamic small- world networks. The network was constructed as a 
2-dimensional Watts-Strogatz model (500 x 500 square lattice with addi- 
tional shortcuts), and the dynamics involved rewiring shortcuts in every 
time step of the epidemic spread. The model of the epidemic is SIR with 
latency time of 3 time steps. The behaviour of the epidemic was checked 
over the range of shortcut probability per underlying bond (j) — ~ 0.5. 
The quantity of interest was percolation threshold for the epidemic spread, 
for which numerical results were checked against an approximate analyti- 
cal model. We find a significant lowering of percolation thresholds for the 
dynamic network in the parameter range given. The result shows that the 
behaviour of the epidemic on dynamic network is that of a static small 
world with the number of shortcuts increased by 20.7 ± 1.4%, while the 
overall qualitative behaviour stays the same. We derive corrections to 
the analytical model which account for the effect. For both dynamic and 
static small- world we observe suppression of the average epidemic size de- 
pendence on network size in comparison with finite-size scaling known for 
regular lattice. We also study the effect of dynamics for several rewiring 
rates relative to latency time of the disease. 



1 Introduction 

The epidemic modelling has become a significant and needed branch of com- 
plex systems research, as we have witnessed the recent epidemic threats and 
outbreaks of human diseases (H5N1 and HlNl influenzas 'TOl |8] or severe acute 
respiratory syndrome [SI [2]) or animal (foot-and-mouth disease [5]) and plant 
diseases alike (e.g. Dutch elm disease [13] or rhizomania [13] )• There are two 
crucial characteristics of the epidemic spread that make it complicated to be 
modelled on the one hand, and costly to be prevented in reality on the other: 
firstly, a number of infectious diseases exhibit long-range transmissions of varied 
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nature, and secondly, the contact network of individuals affected by the disease 
may change in time as the epidemic spreads. These features make epidemiolog- 
ical models a part of larger studies of dynamics on complex networks, but also 
dynamics of complex networks. 

Research findings of the epidemic spread on dynamic networks include its 
behaviour on adaptive networks, where the susceptible are able to avoid con- 
tact with the infected [5]), however a coupling between the epidemic and the 
network dynamics does not necessarily exist. For instance, in [I], spread of 
the aforementioned plant diseases is modelled by vectors performing random 
walk on the network, thus infecting individuals on their paths; Saramaki and 
Kaski fT2] utilise SIR (Susceptible-Infectious- Removed) mechanism on a dynam- 
ically changing small- world contact network, although mainly time development 
of the epidemic is of their interest. Likewise, in [TS] (where focus is on the av- 
erage epidemic size in time) nodes of the contact network can swap their edges 
at a given rate, preserving the degree distribution. It is also worth to note 
[TT] . where disease spread was simulated on a weighted contact network pro- 
duced from real day-to-day encounters (as weights represent the frequency of 
encounters, the dynamics has been in a sense projected onto static weighted 
network) . 

While dynamic network models have been applied in the recent research, 
it seems that we lack comparative study on how the dynamics of the network 
influences the process that takes place on it. The aim of this paper is to find 
and quantify this effect for SIR epidemic spread on static and dynamic small- 
world networks. Based on known analytical calculations for static small- world 
network [7] we derive corrections accounting for the dynamics of the network, 
and check the results against numerical agent-based simulations. 

2 Model 
2.1 Network 

We adopt Watts-Strogatz model of a small-world network 16 : first we take 
a 2-dimensional square lattice with N = L"^ nodes and 2N undirected edges. 
To avoid some finite-size effects we impose periodic boundary conditions for the 
lattice (i.e. we get a torus) . Then, we add a number of undirected edges between 
random pairs of nodes. The number of additional edges ('shortcuts') is set as 
2(j)N, hence (j) is shortcut probability per underlying bond. Network with (j) — 
is just a regular lattice. For nonzero 4> we call the network a static small-world. 

The third type of network is a dynamic small-world. One can construct 
it by randomly distributing shortcuts in every time step of simulation. Here, 
we choose 2(I)N nodes randomly, and keep them fixed for the whole run of the 
epidemic. In every time step we randomly launch shortcuts anchored in these 
nodes, which means the dynamics consists in rewiring one end of these shortcuts. 
For the sake of simplicity we allow for multiple shortcuts being incident with 
the same node, for shortcuts leading to nearest neighbours, and for loops being 
formed. The construction of the source nodes launching shortcuts allows for an 
easier interpretation of the network: the fixed nodes could correspond to centres 
of activity that can be identified as in the real world networks. 
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(a) Regular 2D grid with periodic boundary conditions (torus). |(b)| 
Watts- Strogatz 2D small-world network: 2D grid with shortcuts added to it. 
(c) Dynamic small-world: all the long-range links connected to a set of source 
nodes randomly rewire in time. 



2.2 Epidemic 

The SIR (Susceptible-Infectious-Removed) model is adopted, where the disease 
is transmitted along the edges of the network in discrete time steps. The prob- 
ability p of infecting a susceptible node by an infectious neighbour during one 
time step is set equal for short- and long-range links, both static and dynamic. 
The latency time / of the disease is measured in discrete time units (we take 
/ = 3,4). Thus, an infectious node can transmit disease to susceptible nodes 
with probability p every turn for the period of / turns, and after that time it is 
removed, i.e. it cannot infect nor be reinfected. Every simulation starts with 
only one initially infecting node, all others being susceptible, and it ends when 
no node in the infectious state is left. Sample snapshots of the epidemic time 
development are presented in Figl5] 

Grassberger |3] related the probability of infection to the probability T in 
bond percolation through T = J2t=iPi^ ~ vf^^ = 1 — (1 — p)', where T is the 
so called transmissibility (it is the total probability of a node infecting one of its 
neighbours during the whole latency time) . In the case of 2-dimensional square 
lattice the bond percolation threshold is Tc = 0.5. 



3 Numerical data 

3.1 Parameters of simulations 

The linear lattice size used for most calculations is L = ^/N ~ 500. In Section 
Owe take sizes L = 50, 63, 79, 100, 126, 158, 199, 251, 315, 397, 500. The disease 
latency is set to / = 3 (for faster simulations reported in Sec. [S]) or Z = 4 (in 
Sec. 15.31 in order to get larger set of dynamic rates). The range of probability 
p scanned is p — 0.05 — 0.22 (depending on (f>) with resolution of 1/1024, which 
translates to around T — 0.15 — 0.5. For every p and (j) the epidemic is run 
1024 times with random distributions of shortcuts each time. The fraction of 
shortcuts is (p — — 0.5, with steps of 0.025. The simulations are performed for 
both static and dynamic small-world network. 
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Figure 2: (Colour online) Snapshots of the epidemic spread slightly above 
percolation threshold. L ~ 512, the number of shortcuts is 10 (which gives 
(p = 2 ■ 10~^). t gives the epidemic's time steps. The snapshots for t = 364, 694 
show a dynamic infection (the two joined blue lines appear) . 

3.2 Calculating percolation threshold 

In the study of the epidemic spread on networks, we stick to the percolation 
theory as a reference point. In the theory, a percolation threshold would be the 
value of p that generates an epidemic cluster spanning between the boundaries 
of the whole system. Otherwise, it is possible to define percolation as the point 
at which a cluster of macroscopic size forms (i.e. it occupies a finite fraction of 
the system ioi N ^ oo). We employ the latter to define percolation threshold 
(numerically) as the point at which the average epidemic's size divided by N rises 
above a certain value (here, set to 0.00115). The average is taken over a number 
of reruns for different shortcut drawings. As we can perform simulations only 
for finite sizes, we take the results for a relatively large network of ^/N — 500. 

The choice of the threshold value is taken so as to calibrate the results for 
the static network to the previously confirmed analytical result. We take as the 
theoretical model , where the generating function and series expansion meth- 
ods were used to find the approximate position of bond percolation transition 
in 2D small-world network, which corresponds to the epidemic spread on what 
we refer to as static small-world. 

4 Theoretical analysis 

We can account for the change between static and dynamic networks analytically 
using the model known for static small- world network [7] . As the original theory 
has no time variable, it would be a hard task to introduce dynamics explicitly. 
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The solution, however, is astonishingly simple. One can estimate the average 
number of nodes infected through shortcuts during latency time I: 



Nstat^cf^statN ■T = 4>statN ■Y.Pi^-pf'^ : (1) 

t=l 

i.e. the number of shortcuts in the static network multiplied by the total prob- 
ability of infecting a neighbouring node (this probability is the same for both 
regular links and shortcuts). The analogous expression for the dynamic network 
is found easily 

Ndyn = ^dynN ■ ^ I (^\\^ - pY'" - ^dynN ■ Ip , (2) 

where the sum is an average number of infections transmitted by a single source 
of dynamic shortcuts for a given latency time. It comes from the fact that a 
dynamic shortcut can pass infection several times (the factor p*), while in the 
static case a node could infect only once (since nodes cannot be reinfected in the 
SIR model). This expression predicts lowering percolation thresholds, although 
numerical values of the shift are considerably smaller than the ones obtained 
from simulations. 

Figures 3(a)||3(c) explain why the above expression is not yet correct: it is 



derived only for the source nodes passing the disease on, while it disregards the 
fact that the node may itself become infected via long-range link. Since on the 
static network there is no difference between shortcuts' source and target nodes. 



we can attach the factor (l)N/2 to both infection graphs presented in Fig 3(a) 
For dynamic network, the graphs in Fig 3(b)|3(c)" for infecting a source node 



through a regular link and through a dynamic link give different counts of how 
many shortcuts were used. The former was given in Eq|2] as Ip, and the latter 
actually utilises the same formula, but with the substitution I ^ I + 1. In total, 
we get 

Ndyn = 4>dynN/2 ■ Ip + (t>dynN/2 ■ {I + I) p . (3) 

We assume that Ndyn = Nstat if the epidemic on both networks has the 
same percolation threshold. Thus, we can obtain the ratio of the two shortcut 
densities 

Kp,O^W**..^£^^^^i^ , (4) 



l-(l-p) 



where p is the probability of infection in one time step and I latency time of a 
disease. Now, we can calculate Tc{r(j))) numerically, just as we do it with the 
fitted Tc[(l -I- «)(/>)] in Fig|31 The ratio in Eq|3]was used to plot the lower solid 
line in FigH) 



5 Results 

5.1 Shift of percolation thresholds 

In Figure|4]we plot numerical and theoretical values of percolation thresholds Tc 
for both static and dynamic small-worlds. The resulting Tc((/)) data points for 
static small-world network agree with the analytical approximation 7 , which 
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Figure 3: (a) Infections through static shortcuts are symmetric. |(b)| Infection 



of the dynamic shortcuts' source through regular lattice, 
dynamic shortcuts' source through a shortcut. 



(c) Infection of the 



confirms the validity of calibration procedure. As the lower dataset marks the 
effect of network dynamics, the difference between the two networks proves 
to be systematic and significant. The dashed line is a fit Tc[{l + v)(p] of the 
analytical model for the static network, where the fitted parameter v may be 
interpreted as a virtual percentage of additional shortcuts needed to obtain the 
dynamic network percolation thresholds. It follows from the fit that percolation 
thresholds for dynamic network are lower as if the shortcut density were {l + v)<j) 
(where v = 0.207 ± 0.014 is the fitted parameter). Nonetheless, qualitatively 
the epidemic on dynamic small world behaves in the same way as on the static 
one for the given range of parameters (0 — 0.5 corresponds to every node in the 
network having on average two additional links). 

The analytical correction slightly exceeds the values of simulation data points, 
but the overall agreement is satisfactory. The difference between the analytical 
solution and the observed behaviour does not exceed the shift between static 
and dynamic networks obtained from simulations. The discrepancy might be 
due to the method of calculating percolation thresholds from numerical data or 
due to the approximate nature of the correction. 



5.2 Suppression of finite-size scafing 

The primary motivation of checking finite-size scaling for the system was to 
utilise it to determine the percolation thresholds very accurately (as the shift of 
thresholds observed in Fig|3]is relatively small), and to arrive at threshold value 
for infinite system size. Yet, it is worth noting at this point that the knowledge 
of thresholds for infinite system sizes would not usually be appropriate for eval- 
uation of risks in the real epidemic, given the sizes of some real networks. To 
study the size of finite-size effects is thus vital on its own right. 



In Figure 5(b) the convergence of the average epidemic size to the threshold 
behaviour can be observed, and the significant dependence on system size ranges 
up to the epidemic size of around 0.57V and interval of transmissibility of length 



around 0.08 (the numbers are very rough estimates). As presented in Fig 5(a)[ 



one may check that sections of the plot for a given average epidemic size obey 
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Figure 4: (Colour online) Circles dataset - the static small- world. Squares 
dataset - the dynamic network. The solid blue (upper) line is the analytic 
approximation [7] for Tc(0) and the dashed line gives Tc[{l +«)(/))], with the fit 
parameter v — 0.207± 0.014. The solid purple (lower) line represents theoretical 
approximation from Sec. 15.31 Error bars are of the size of the plot markers, 
unless visible. 
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scaling of the form 

Tn^T^- N-'/'' = Too - i"'/" , (5) 

where Tn are the values of transmissibility for a given system size N and a set 
section position, and Too is the percolation threshold for infinite system size. For 
regular lattice Too is fitted correctly for various section positions as 0.500±0.005 
(the error may vary for different sections, but does not exceed the given value). 

It appears that the dependence on system size for small-world networks 
(both static and dynamic) is dissimilar to the one of regular lattices, as can be 
seen in Fig l5.2l {(/> — 0.05). It is suppressed to smaller values of the average 
epidemic size. For the shortcuts density (p = 0.5 the dependence on system size 
is already visible only below the epidemic size of 0.03. Because the dependence 
of the epidemic size on size of the system becomes of the order of magnitude 
of statistical fluctuations (the quality of the data can already be seen in the 
Fig l5.2p . any attempts to utilise finite-size scaling for determining percolation 
threshold are not viable. Indeed, the errors do not allow us to check if the same 
form of finite-size dependence as in Eq|5] holds. 

5.3 Dependence on the rate of dynamics 

One can generalise the theoretical analysis for various rates of dynamics, given 
the formula in Eql^l To explain this, let us notice that there are two time 
scales in the model: the latency time I of the infection and the duration 1/d 
between consecutive rewirings of dynamic links (both measured in discrete time 
steps of the epidemic spread). As the choice of latency I only rescales the 
total probability of infection T = T(j),l), we can dispose of it, and the crucial 
parameter Id that accounts for the shift of percolation thresholds is defined as 
the number of shortcut movements during latency time. 

Obviously, for a static network we get d = 0, while for all the above analysis 
of dynamic network we have Zc? = 3 (/ = 3 and the rewiring was performed every 
turn, so d = 1). Depending on the interpretation of the model, we could also 
consider d > 1. However, if p is to be the probability of infection during one 
time step it is reasonable that shortcuts rewiring faster than one time step would 
infect with appropriately smaller probability, and there would be no further shift 
of percolation thresholds. 

Since the epidemic spreads with discrete time, which results in sums as in 
Eql21 we are interested in rational numbers d £ [0, l]nQ, particularly of the form 
l/i,i S Z. What we need is N^yn calculated in a similar way to that in Eql^] 
Here, we take I = 4, d = 1, 1/2, . . . , 1/7, and we plot both the numerical and 
theoretical results for <j) — 0.25 in FiglT] Theoretical derivation is to be found 
in the Appendix. The theoretical approach gives slightly exceeding values (the 
scale should be noted), which is the same effect as discussed at the end of Section 

6 Discussion 

We have shown that introducing dynamics of the long-range links in a small- 
world network significantly lowers an epidemic threshold in terms of probability 
of disease transmission, although the overall dependence on number of shortcuts 
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Figure 5: The behaviour of the epidemic outbreak magnitude for various system 
sizes (hnear size vary between L = 50 — 500, left- and rightmost data points in 
|(b)[ respectively), (a) Finite-size scaling Tjv = To^ — L^^^'^ on regular lattice. 
The points correspond to values of T at the level of the epidemic size 0.1. |(b)| 
The extent of size dependence for regular lattice. 
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Figure 6: For dynamic small world size dependence of the epidemic outbreak 
magnitude is suppressed. Inset shows enlarged region around percolation thresh- 
old. 
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Figure 7: Dependence on dynamics for (f> — 0.25, latency I — 4. 

stays the same. Consequently, the risk of an epidemic outbreak is higher than 
in any calculations involving static models. The effect remains secondary to the 
influence that the introduction of additional of shortcuts has on the spread of 
the disease. It should be noted that the shift of percolation thresholds depends 
on the relative measure of dynamics of the network with respect to the process 
on the network (rewiring rate and latency time, respectively). Any accurate 
analytical calculation or simulation should take this quantity as a significant 
parameter, to be estimated for a particular disease and type of the network. 

As in reality we consider only finite-size networks, and real epidemic sizes 
do not usually reach values of the order of even 10% of the system size, the 
information on finite-size effects seemed very much needed. That the epidemic 
outbreak magnitude does not depend on the system size for small-world net- 
works as much as it does for regular lattices means that we should not expect 
the epidemic outbreaks below transmissibility threshold value. Thus, finite-size 
effects seem to become secondary, as well. 

The usefulness of such a model for risk prediction still depends on our knowl- 
edge of the probability of transmission {p or T) of a given disease, which is not 
easy to obtain for diseases spreading outside of well controlled environments like 
hospitals. Relatively good estimates, thanks to the nature of transmission, exist 
for syphilis. Transmissibility of the disease is reviewed in 3, , where authors give 
values ranging from 9.2% to 63% per partner, and decide on 60% as the lower 
boundary for untreated disease. This seems to be well above the epidemic thresh- 
old, irrespective of very different network topology for such diseases. However, 
this also shows that errors on estimates of transmission probabilities exceed the 
effect of threshold shifting studied here. 

Though the 2-dimensional network structure used here may be said to corre- 
spond mainly to that of plantations, it is worth noting its generality: nodes may 
be interpreted as plants, animals or humans, but also on a larger scale as farms. 
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households, or cities and airports; in turn, long-range links could mean wind 
(on farms), disease vectors, occasional human contacts, or airline connections. 
Still, it has some other fairly realistic characteristics: according to [TT], who 
analysed the structure of human social interactions, 'the majority of encounters 
(76.70%; 75.26-78.07) occur with individuals never again encountered by the 
participant during the 14 days of the survey.' This may mean that about 24% 
of the repeated contacts corresponds roughly to our regular underlying lattice 
with z = 4 neighbours for each node, while the 76% correspond to around 3z dy- 
namic contacts distributed over over 14 days. This gives on average w 0.20 for 
simulation with daily time steps, which lies within the parameter range studied 
in this paper. 
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Appendix: Dependence on the rate of dynamics 

Below we present the way to calculate Ndyn for latency periods I — A,b (in the 
simulation we set I = 4, but we need to take into account also the process from 



Fig 3(c) which in a sense increases latency by 1). Let us define 



Ao(p,0 = i-(i-p)' = r(0, i>i 

Ai{p,l) = T{l)[l-T{l-l)] + [l-T{l)]T{l-l) + 
+ 2T{p,l)T{l-l) 

A2{pj) = r(2)[i-r(/-2)] + [i-r(2)]T(z-2) + 

+ 2T{2)T{l-2) 

An{p,l) - 3T{lfT{l~2) + 2T{lf[l-T{l~2)] + 

+ 2T{1)T{1-2)[1-T{1)] + 

+ T{l-2)[l-T{l)f + (6) 

+ 2r(i)[i-T(i)][i-r(/-2)] 

Ai2{p,l) = 3r(l)T(2)r(/-3) + 2T(l)r(2)[l-T(/-3)]-H 
+ T{l)[l-T{2)]T{l-3) + 

+ [i-r(i)]T(2)r(/-3) 

+ r(i)[i-r(2)][i-T(?-3)] + 

+ [l-r(l)]T(2)[l-T(?-3)] + 

+ [i-r(i)][i-T(2)]r(/-3), 

where we substituted r(l) for p on the right-hand sides, and we leave out 
the argument p in T{p, I) to simplify the notation. Those quantities correspond 
to the average number of infections during one latency period depending on 
when the rewiring takes place. One can present those diagrammatically (here 
for Z = 5) as 
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A ob,5) = 

A i(p,5) =.|.... + ....|. = 2.|.... 

A 2b,5) =..|... + ...|.. = 2..|... (7) 

A n(p,5) 

A i2(p,5) =.|..|.. + ..|..|.=2.|..|.. 

where the symbol '|' marks rewirmg, and '•' one epidemic time step during 
latency period. For instance ■\ ■ ■ would correspond to three turns with one 
rewiring, during which either 0, 1 or 2 infections are possible. The derivation 
involves only very easy combinatorics, but for longer latency periods one would 
need to repeat these calculations to obtain more terms and different prefactors. 

Now, one can easily obtain expressions for Ndyn for any l/d G Z. Below we 
give only the general expression for 1/d > I: 

Ndyn = ^^d{[2Ai(/)+A2(/) + (l/d+l-0Ao(0] + 

+ [2Ai{l + l) + 2A2{l + l) + {l/d-l)Ao{l)]} (8) 
where / = 4. The first term in the brackets corresponds to Fig |3(b)| and the 



second to Fig 3(c) For greater numbers of rewiring per turn d, we need to 
consider the terms Aii,Ai2. The result is plotted against simulated data in 
FiglZl 
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